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Abstract. We address in this paper the study of a geometric evolution, cor- 
responding to a curvature which is non-local and singular at the origin. The 
curvature represents the first variation of the energy M P (E) defined in ( |1.I[ |, 
proposed in a recent work [5] as a variant of the standard perimeter penaliza- 
tion for the denoising of nonsmooth curves. 

To deal with such degeneracies, we first give an abstract existence and 
uniqueness result for viscosity solutions of non-local degenerate Hamiltonians, 
satisfying suitable continuity assumption with respect to Kuratowsky conver- 
gence of the level sets. This abstract setting applies to an approximated flow. 
Then, by the method of minimizing movements, we also build an "exact" cur- 
vature flow, and we illustrate some examples, comparing the results with the 
standard mean curvature flow. 



1. Introduction 

In a recent paper [5], the last two authors, together with M. Barchiesi, S. H. Kang 
and T. Le, proposed a variational model for (binary) image denoising, which was 
supposed to preserve small scale details (or small oscillations of the boundary) while 
regularizing the large scales. This model is a variant of the celebrated Mumford- 
Shah functional, where the perimeter term is replaced by the following one 

(1.1) M P {E) = \{x e R d : dist(a;,-B) < p ,dist(x,R d \ E) < p}\ , 

Zp 

defined for any E C R d , where p > acts as a scale selection parameter. Here and 
throughout the paper, given A C R d measurable, we denote by \A\ its Lebesgue 
measure. Notice that the energy is finite if and only if dE is compact. The idea 
behind such a variant is that fluctuations of dE at lengths much smaller than p > 
will have very little influence on the energy; on the other hand, it behaves as the 
standard perimeter on much larger smooth boundaries, and in a more complicated 
non-local way on sets with fine microstructures on scales of order p. 
The sort of denoising which is obtained in [5] is shown in Fig. [Tj where small 
oscillations (here the stripes of the fingerprints) are almost untouched, while the 
noise has mostly been removed. 

In this paper, we try to investigate some mathematical analysis aspects of this 
model. More precisely, we want to study the geometric evolution of curves and 
shapes by the gradient flow of the functional proposed by these authors. 

To this purpose, we first extend our energy to L 1 functions, and express it in 
terms of a function depending on the oscillation of u on balls of radius p, following 
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the approach in [16] . With this point of view, it turns out that (1.1 ) is the restriction 
to characteristics functions of a convex, l.s.c. functional, satisfying a suitable "coarea 
formula" . Then, we introduce the "curvature" as the first variation of this functional 
with respect to inner variations of the sets. This curvature is not continuous and it 
is not well defined for all smooth sets. Therefore, in (2.6) we introduce a smoother 
version M.f of (1.1), which roughly speaking consists in averaging M. r over r, for 
r varying in a neighborhood of p. The corresponding curvature is now well defined 
on smooth sets. 

After this preliminary analysis to define a proper notion of curvature, we study 
the corresponding geometric flow. Using a level set approach and working in the 
framework of viscosity solutions, we define a mean curvature flow equation, which is 
both non-local and singular. Indeed, our Hamiltonian F(x, Du, D 2 u, K) depends in 
a non-local way on the level set K, and behaves like a power (d— 1) of the curvature 
tensor of dK for vanishing sets, being thus singular in dimension d > 3 (see (3.27)). 
To deal with such degeneracy we combine the approach by Slepcev [26] to non-local 
Hamiltonians with the approach by Ishii and Souganidis [53] and Goto [2D] to de- 
generate Hamiltonians. However, the approach in [55] is based on the assumption 
that the Hamiltonian is continuous with respect to all its variables, in particular 
with respect to L 1 convergence of the sets. This is not the case of our Hamiltonian 
(and of any reasonable regularization of it). Therefore, we build up a variant of 
the approach in [5B] that works for a general class of Hamiltonians satisfying suit- 
able continuity properties with respect to the Kuratowski convergence instead of 
L 1 convergence of sets. We adapt the notion in [23] of viscosity solutions for sin- 
gular Hamiltonians to our non local setting, and we show a corresponding result of 
existence and uniqueness. This result will apply to a general class of Hamiltonians, 
which does not include the Hamiltonian corresponding to our non-local curvature 
flow, but only a suitable continuous approximation of it. 

Finally, we study the minimizing movements corresponding to the energy . 
We introduce an implicit time-discretization of the motion, and we show that it 
converges, up to a subsequence, to a solution of the level set equation in the vis- 
cosity sense. In this way we recover an existence result for viscosity solutions also 
for the exact Hamiltonian corresponding to the first variation of , yet without 
uniqueness. We mention that in a recent paper of Caffarelli and Souganidis [9], a 
similar strategy (based, this time, on a diffusion/thresholding time-discrete scheme) 
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has been successfully implemented to build up a non-local curvature flows associ- 
ated to fractional diffusions. Our time-discretization approach can be numerically 
implemented, following the approach in |16j : we show eventually in Section [5] a few 
examples which are compared with the classical mean curvature flow, and seem to 
confirm a slower smoothing of oscillatory boundaries. 

We mention the existence of a few interesting alternative approaches to non-local 
evolutions. The recent papers of [5TJ [5] provide a point of view slightly different 
from ours, and address different kinds of evolutions. In particular, [8j also deals with 
non-monotone evolutions, such as the one describing the motion of dislocation lines 
in crystals (see also [3] |H1 [33] ) . Another approach is described in the papers 
of Cardaliaguet [TUl [TT] , Cardaliaguet and Rouy [TS] , Cardaliaguet and Ley [121 Q3] ■ 
There, appropriate definition for evolving tubes are proposed and the convergence 
of a time-discrete scheme (of the same kind as ours) is addressed in [14J . Moreover, 
except in the preliminary work |10j , the authors of this series of papers have taken 
care to never need to evaluate the velocity on arguments which are not "natural" 
(such as smooth level sets and their normal or second fundamental form) , contrarily 
to what is needed in our proof of uniqueness (as in |26j). Unfortunately, extending 
their work to our approach raises complicated technical issues, since in particular 
our velocity does not have the required continuity properties, and our minimizers 
have unknown regularity. This is an interesting direction for future research, but we 
also believe it is useful to develop the level-set approach in the non-local geometric 
setting. 

To summarize, the first goal of this paper is to investigate the geometric flow 
corresponding to a non-local variant of the perimeter introduced in [S] , in connection 
with image denoising. We have developed a viscosity approach to non-local singular 
Hamiltonians, combining many ideas from [26], [23] and [20]. Through the viscosity 
approach we have obtained existence and uniqueness for a suitable regularization of 
our Hamiltonian, while a minimizing movements approach yields a solution for the 
original Hamiltonian. The abstract approach is itself interesting and stimulating: 
a complete picture at the moment is still missing, and this paper represents a first 
attempt to study singular non-local Hamiltonians, not even continuous with respect 
to L 1 convergence, but only with respect to Kuratowski convergence of level sets. 
The combination of the minimizing movements variational method with viscosity 
techniques seems to be a promising approach to complete the picture. To our 
knowledge, up to now this kind of study has been carried out only in |14j . |18) . 
and a few papers by the first author and co-workers. We hope that (borrowing in 
particular from [T3]) we will be able to extend these ideas to other motions, and 
also, understand how to make the proof of the comparison result less dependent on 
the extension of the Hamiltonian out of its natural domain of definition. 



2. The energy functionals 



2.1. The p-neighborhood. As mentioned, we focus on the study of ( |1.1[ ). It is 
well-known that, under mild regularity assumption on E (see for instance [3]) we 
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have 



UmMJE) = H (dE) = Per(E), 

p^>0 



where Per(E) is the standard perimeter of E. It is also very easy to show that A4 P 
T-converges to the standard perimeter [17]. An issue with definition (1.1 ) is that it 
depends on the choice of the representative within the Lebesgue equivalence class 
of the set E. For this reason, one may introduce the following variant: 

(2-1) £ P (E) = — I osc B{XiP) {x E )dx 

where osc^(m) denotes the essential oscillation of the measurable function u over a 
measurable set A, defined by 

oscyi(u) = ess sup u — essinfw. 

A A 

One checks that £ p {E) coincides with the measure of the p-neighborhood of the 
essential boundary of E. Moreover, 

£p(E) = m{{M p (E') : \EAE'\ = 0}, 

where EAE' denotes the symmetric difference (E\E')U (E' \ E). Finally, £ P {E) = 
£ p (E c ) and it is finite if and only if either E or E c is (essentially) bounded (where 
E c := R d \ E). An advantage of Definition (2.1 ) is that can be easily generalized to 
a measurable function u € Lj oc (R d ). By a slight abuse of notation, we still denote 
£ P {u) the functional: 

(2.2) £ p {u) = o / osc B(x , p) (u)dx. 

One can check that this energy is one-homogeneous, convex and therefore sub- 
additive. Moreover, it is lower semicontinuous with respect to weak* convergence 
in Lf£ c , and satisfies the following generalized coarea formula 

/oo 
£ p {{u > s})ds. 
-OO 

This follows from the fact that for any A and u we have 

/oo 
osc A {X{u>s})ds. 
-oo 

Moreover, one easily deduces that, given two measurable sets A,Bc K d , 

(2.4) £ P (AUB) + £ P {AC\B) < £ P {A) + £ p {B) . 

Indeed, observing that xaub + XAnB = XA + Xb, by coarea formula and in view of 
the subadditivity of £ p , we have 

/oo 
£ P ({XAuB+XAnB > s})ds 
-OO 

(2-5) = £ p (xa + Xb) < £ p (xa) + £ p (xb) = S„(A) + £ p (B). 
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2.2. The continuous energy functional. Fix p > and 6 € (0, po). We consider 
now a Lipschitz function / : R — > R + which is even, with supp/ = [— po,Po], 
constant in [— <5, 5] and nonincreasing in ]R + . 



We then introduce the following variant of ( 1.1 1 



(2.6) M f (E) = / f(d E (x))dx 



where d E is the signed distance to dE (negative inside E and positive outside). 
Notice that M.* (E) is finite if and only if dE is compact, i.e., E or its complement 
is bounded. We now show that 



fPO 

M f (E) = / (-2sf(s))M s {E)ds. 
Jo 



ii 



Since M^(E) — (E c ) we can assume that E c is bounded. Then, thanks to the 

co-area formula we have 

r rpo 
M f (E) = / }{d E {x))\Dd E {x)\dx = / f{s)H d - l (d{d E > s})ds 

Pa 

-f'( s )\{s <d E < p }\ds 

Pa 

Pa 

-/'0)(|{- s <d E < poll - \{s <d E < p }\) ds 

Pa 

{-2sf'(s))M s (E) ds. 

o 

Thanks to this, we can introduce the following variant of , defined on Borel sets 
and which depends only on the Lebesgue equivalence class 

fPO 

(2.7) £ f (E) = / (-2sf'(s))£ s (E)ds, 

Jo 

As before, we consider the convex extension of , defined for all functions u € 
Lj oc (R d ) by 

rPo r rpo 

(2.8) £ f (u) = / (-2sf{s))£ s (u)ds= / / (-2s/'(s))osc B(x , s) (u) dsdx. 

Jo JnJo 

By construction, £f is a convex, lower semicontinuous energy which satisfies the 
generalized coarea formula 

/•oo 

(2.9) £ f (u) = / £ f ({u > s})ds. 



Clearly, (2.4) is still true also for £■? . 



2.3. The non-local curvature. Let E C K d be a smooth set with compact bound- 
ary. We denote by v E (x) the outer normal unit vector to dE at x. The non-local 
curvature k p is formally defined as the first variation of the energy £ p in (2.1 ). Set 



(2.10) k p (E,x) = k+(E,x) + k-(E,x), 



6 



A. CHAMBOLLE, M. MORINI, AND M. PONSIGLIONE 



where 



(2.11) k+(E,x) 



2p 



det(I + pVve{x)) if dist(a; + pve{x), E) = p, 



(2.12) k~(E,x) 



2p 



det(I - pVi/ E (x)) 



otherwise, 

if dist(x — pve{x), E c ) = p , 
otherwise. 



Let us decompose dE into three sets: dE = A+ U U Af^~, where 

:= {x E dE : there exists t > p : dist(x + tv E {x),E) = t}, 
B+ := {x e dE : dist(sc + pv E {x),E) < p}. 
Af+ = dE\(A+UB+). 

Analogously, we define^4~ , B~ , Af~ with u E (x) replaced by —v E {x), and we set 
AA,:=AA+UAA-. 

Lemma 2.1. Let E C W 1 be a set of class C 2 with compact boundary, such that 
H^iAfp) = 0. Then, for every ip € C 2 (dE;R) we have 
d 



(2.13) 



K p (pd'H 



d-i 



where <& e is a diffeomorphism such that $(x) = x + E(p(x)ve{x) for x € dE. 



Remark 2.2. Notice that the assumption of Lemma |2.1| holds true for generic 
smooth sets. More precisely, given a smooth set E, then for almost all positive p 
one has T-L d ~ 1 {N p ) = 0. Moreover, such assumption is crucial. Indeed, let E be a 
rectangle of sides 2 and 4, respectively, and set p = 1. In this case, a curvature k p 
satisfying (2.13) is not well defined. Indeed, one readily sees that such curvature k p 
should depend in a non-local way on (p. More precisely, let ip = (pi + ip2, where tpi 
are defined in a small neighborhood of the middle points Pi of the large sides of 
E, and assume that tpi have constant sign. Then, if + ^2(^2) > 0, then (2.13) 

holds true, while if yi(pi) +^2(^2) < 0, then (2.13) holds true with n p replaced by 
Kj. In particular, the "curvature" at p\ depends on the value of ip at P2- 



£+{F\E) 



Proof of Lemma \2.1\ We will show that 

(2.14) ^£f(^(E)\E) u __ o 

where, for every set F 
1 

2p ./ 



n±pdU c 



°8c B (x, p )(xf) dx and £ p (F\E) 



1 

2~i> 



osc B{x . p) (xF)dx. 



(2.15) 



In order to prove (|2.13|), we will focus on the identity 

K+<pdH d -\ 



de 



i)E 



the variation of £ p being analogous. The proof is divided in three steps: first, we 
prove that (2.15) holds if the support of (p is contained in A+. Then, we show that 
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the variation vanishes on Bj. Finally, by a localization argument, recalling also 
H d ~ 1 (N+) = 0, we deduce the validity of (|2~15l. 



Step 1. In this Step we assume that supp(ip) C At, and then prove (2.15). For 
every x € let y s (x) := $ e (x) = x + £y(cc)^£;(:E), and set 

(2.16) E+ e := {y s (x) + *i/* e(JS) (ife(aO), x£dE,t£ (-p,p)} \ E, 

so that, for e small enough 

f+(* 6 (^)|£) = \E+ e \. 

Set now 

£+ e := {x + tv E (x),x €dE,te (0, p + e(p(x))} \ E. 

By construction, since dE is of class C 2 , one can see that \E+ e AE+ £ \ — o(e). 
Finally, we have 

£+(* e (E)\E)-£+(E\E) |£+ £ |-|<ol = \E+ E \-\E+ fi \ | 



e 



dx / |det(/ + (/o + f)Vi/ E (x))|dt + o(l). 



SB JO 



For e — > we recover (2.15) 



Step 2. In this step we show that the curvature k+ vanishes on This amounts 
to show that, if tp has support in then 

This is readily seen, since by definition of £?+ we have that, for e small enough, 
E+ e = E+ , so that 

E+(* e {E)\E) = \E+J = \E+ \=£+(E\E). 



Step 3. In this step we conclude the proof by standard localization arguments. 
Given S > 0, we can always write p = pi + p>2 + P>3 where supp(<^i) C A+, 
supp((,C2) C ■H d_1 (supp(( / 3 3 )) < S and \pi\ < \<p\. Notice that 

£+(<Z> e (E)\E) = \E+ i+V2+ip3 J = \E+ +lp2 J+r, 

where |r| < CeS. Moreover, for e small enough 

Therefore, using Step 1, and since n^P2 = by Step 2, we conclude 



(2.17) 



lim sup 

£->0 



£+($ s (E)\E)-£+(E\E) 



< 



lim 

£-S-0 



£ JdE 

We conclude by the arbitrariness of 6. 



K+p 1 dU d - 1 



C'S = C'S. 



□ 



A. CHAMBOLLE, M. MORINI, AND M. PONSIGLIONE 



Now we introduce the non-local curvature Kf associated to the energy (2.7): 



(2.18) 
where 
(2.19) 



Kf(E, x) — Kf(E, x) + Kj (E, x), 
nf(E,x)= f\-2sf'{s))4{E,x)ds. 



The following result is a direct consequence of (2.7), Lemma 2.1 and Remark 2.2 



Theorem 2.3. Let E C 

ip G C 2 (dE;R) we have 



be an open bounded set of class C 2 . Then, for every 



(2.20) 



d 



d^ £ ^ E K- 



K f {E,x)(f(x)dn d - 1 (x) 



i)E 



where <I> e is a diffeomorphism such that $(x) = x + elp(x)ve{x) for x € dE. 

3. Viscosity solutions of the non-local level-set equation 

In this section we introduce the level set formulation of the geometric evolution 
problem 



(3.1) 



V = K f , 



where V represents the normal velocity of the boundary of the evolving sets t h-> E t , 
and we give a proper notion of viscosity solution. Then, we develop an abstract 
setting where we provide existence and uniqueness results, for a suitable class of 
Hamiltonians. Unfortunately, as already said in the introduction, this theory applies 
only to a suitable regularization of the curvature Kf. However, we will also provide 



later on (Section Bl) an existence result, without uniqueness, for equation (3.1). 



3.1. The non-local evolution. Here we introduce the level set formulation of the 



geometric evolution problem (3.1). To this aim, following the level set approach, 
we identify E t with the superlevel set {u > 0} of a function u : K x R d ^ K, and 
study the corresponding degenerate parabolic equation in u. Let Mf y x ^ denote the 
class of d x d symmetric matrices, and C(M. d ) the class of closed subsets of R d . We 



introduce the Hamiltonian F, 



f 



tidxd 



(3.2) 
where 



„ vm x C{R d ) i — ^ R defined by 

fPO 

F f (x,p,X,K):= / (-2sf(s))F s (x,p,X,K)ds, 
Jo 



(3.3) F s (x,p,X,K) = F+(x,p,X,K) + F-(x,p,X,K) 

and the functions F c + and Fr arc defined as follows: 



F+(x,p,X,K) 



bl 

2s 




det 



I 



]p\ VpXVp . 



if 



dist(a; — sp, K) > s, 
otherwise, 
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F-(x,p,X,K) = 



2s 



if 



dist(a; + sp, K c ) > s, 



\P\' 

otherwise. 

Here := {I — p ® p), where, for p / 0, p = p/IpIi an d f° r a symmetric 
matrix, is the matrix with all eigenvalues replaced with their positive part (in 
particular, det[X] + = for any X which is not positive definite). 

Remark 3.1. If u is a smooth function and u{x) is not a critical level, by Theorem 
|2.3| we easily deduce that 



(3.4) F f {x,Du(x),D 2 u,{y e R d : u(y) > u(x)}) = \Du(x)\k { {{u > u(x)}, x). 

In this identity we use in particular that, if p = Vit, X = D 2 u, K = {u > u(x)}, 
r 1 + 

then det I — ArVpXVp = means that there is a direction along which the 
curvature is larger than 1/s, so that k^{K,x) = 0. 

The level set approach consists in solving the following parabolic Cauchy problem 

{u t (x,t) + F f (x,Du(x,t),D 2 u(x 7 t),{y: u(y,t) >u(x,t)}) =0 
for t > 0,x e R d , 
u(0, •) = uoi 

in the viscosity sense. The definition of a viscosity solution for such a non-local 
Hamiltonian will be introduced in the next subsection. We will prove an existence 
and uniqueness result in this setting, which will be applied to a smoothed variant 
otF f . 

3.2. The abstract setting. We introduce here a notion of viscosity solutions for 



problems such as (3.5 1. The issues are of course that the Hamiltonian is non- 
local, but also that it is singular in p — (at least in dimension d > 3), in the 
sense that it grows as the set vanishes as a power (d — 1, in dimension d) of the 
curvature tensor. For this reason, we have to adapt both the setting of Slepcev 
[26j for non-local evolutions (notice however that we will consider weaker continuity 
assumptions with respect to the set variable), and the one of Ishii and Souganidis 
[2"3"] (see also Goto [5D]) for singular Hamiltonians. 

We will first list the properties which our Hamiltonians need to satisfy in order 
to show an existence and uniqueness result, and then introduce the appropriate 
definition of a viscosity solution (which is almost standard). Let A,(Mr) denote 
the family of open sets in M. d , and C(M rf ) the family of closed sets. We consider 
Hamiltonians F : R d x R d \ {0} x M^f x {C(R d ) U A{M. d )} i-> K satisfying the 
following properties: 

i) Translational invariance: Fix + r, p, X, E + r) = F(x,p,X,E) for every 
r e R d ; 

ii) Degenerate ellipticity: F{x,p,X,E) > F(x,p,Y,E) if X < Y; 

iii) Monotonicity in the set variable: F(x,p, X, E) > F(x,p, X, G) if E C G; 

iv) Geometric property: F(x, Xp, \X+[ip®p, E) = XF(x,p, X, E) for all A > 0, 
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v) Continuity: F is continuous with respect to its first variable, moreover, the 
following properties hold: 

v.f) If x n —> x, p n —> p 0, X n X and {K n } C C(]R d ) is a sequence 
converging to K in the Kuratowski sense, then 

F(x,p,X,K) < liminf F(x n ,p n ,X n ,K n ). 

n 

v.2) If x n — > x, p n — > p ^ 0, X n — > X and {A n } C .4(M d ) is a sequence 
such that A n converges to A c in the Kuratowski sense, then 

F(x,p,X,A) > limsup F(x n ,p n , X n , A n ). 

n 

vi) There exists a continuous function c : (0, +oo) i— > (0, +oo) such that, for all 
x eR d , p EM. d \ {0}, E e C(R d ) U A(R d ) we have 

(3-6) - c(\p\) < F(x,p,±I,E) < c(\p\). 

Following [23], we introduce the family T of functions / G C 2 ([0, oo)) such that 
/(0) = /'(0) = /"(0) = 0, and such that f"(r) > for all r > which satisfy 

(3.7) hm/Mc(b|)=0. 

We refer to [UJ p. 229] for the proof that the family T is not empty. 

Let T > be fixed. As a slight variant to [23], we introduce the following 
definition. 

Definition 3.2. We will say that if G C°(R d x (0, T)) is admissible at the point 
I = (x,t) £/ it is of class C 2 in a neighborhood of z and, in case Dip(z) = 0, the 
following holds: there exists f £ T andcu G C ([0,oo)) satisfying lim r _ v o tL, ( r )/ r = 0j 
such that 

\rtx,t)~ V (z)~<p t (z)(t~i)\ < /(|a;_f|)+a;(|t-t|) 
for all (x,t) in a neighborhood of z. 

Given a function uq, which is uniformly continuous in M. d , we want to solve 

(3.8) ih(x, t) + F(x, Du(x, t), D 2 u{x, t), {y : u(y, t) > u(x, t)}) = 

for (x,t) G R d x (0,T), 

subject to the initial condition it(0, •) = uq. We introduce the following definition 
of a viscosity sub/supersolution, inspired from both frameworks of [23] and [26] . 

Definition 3.3. An upper semicontinuous function u : R d x [0,T) — > K is a viscosity 
subsolution of if for all z := (x, t) G R d x (0, T) and all ip G C°(R d x (0, T)) 

such that u — ip has a maximum at z and if is admissible at z we have 

[ ip t {z) + F (x, Dip(z),D 2 <p(z), {y : cp(y, t) > cp{z)}) < if D<p(z) + 0, 

(3.9) <^ 

I ft{z) < otherwise. 



A lower semicontinuous function is a viscosity supersolution of (3.8) if for all z G 
M. d x (0,T) and all ip G C°(Mr x (0, T)) such that u — p has a minimum at z and 



A NON-LOCAL MEAN CURVATURE FLOW 



11 



tp is admissible at z we have 

<p t (z) + F (x,Dtp(z),D 2 tp(z), {y ■ <PM > <p(z)}) > ifDtp(z) ? 0, 
¥t{z) > otherwise. 



(3.10) 



Finally, a function u is a viscosity solution of (3.8| if its upper semicontinuous 



envelope is a subsolution and its lower semicontinuous envelope is a supersolution 



of (3.8) 



As it is standard in the theory of viscosity solutions, the maximum in the 
definition of subsolutions can be assumed to be strict, while the test functions tp 
can be assumed to be coercive (and similarly for supersolutions) . For the reader's 
convenience, we show that this is the case also in our non-local setting. Assume 
for instance that u is a subsolution, u — tp has a maximum at some (x,t), with tp 
admissible at (x,t). If Dtp(x,t) ^ we replace tp with 

<Pe(y, s ) : = <p(v, s ) + e \y - x \ 2 + 1* - s l 2 - 



Then the maximum of u — tp e at (x,t) is strict, and we recover the inequality (3.9) 
for ip by letting e — > and using the semicontinuity of F, observing that the sets 
{tp(y,t) + s\y — x\ 2 > ip(x,t)} converge to {ip(y,t) > (p(x,t)} in the Kuratowski 
sense. We use then the semicontinuity property v.l) to conclude. 

If Dip(x, t) — 0, we choose / G T as in Definition |3.2| and replace ip by 

<p{y, s ) ■= <p(v, s ) + f(y -x) + \t- s\ 2 . 

We still have D(p{x, t) = 0, dp is admissible at (a;, t), (f>t{x, t) = <pt(x, t) and now the 
maximum of u — dp is strict. 

Notice that our definition of supersolutions and subsolutions is formally different 
than the one given in [55], that involves the superlevel sets of u instead of tp. Indeed, 
in the case of a subsolution we can assume that the test function tp is such that 
U < tp, and u(x,t) = tp(x,t). Then, {y : u{y,t) > u(x,t)} C {y : tp(y,t) > tp(x,t)} 
so that 

F (x, Dtp(x, t),D 2 tp(x, t),{y: u(y, t) > u(x, t)}) 

> F (x, Dtp(x, t),D 2 tp(x, t), {y : tp{y, t) > tp(x, i)» . 

Therefore, our definition seems actually weaker. The following Lemma shows that, 
in fact, it is equivalent. 



Lemma 3.4. Let u be a viscosity subsolution of (3.8). Then, for all (x,t) in 
R d x (0,T) and all tp € C°(R d x (0,T)) admissible at (x,t), with Dtp(x,t) ^ 0, and 
such that u — tp has a maximum at (x,t) we have 

(3.11) tp t {x,i) + F (x,Dtp{x,t) 1 D 2 tp(x,t),{y : u(y,t) > u{x,t)}) <0. 
A similar statement holds for supersolutions. 

Proof. We can assume that the test function tp is such that u < tp and u(x, t) = 
tp{x,t). Consider a decreasing sequence tp n of functions which are smooth and 
such that inf„ ip n — u, ip n > u + l/n. Such a sequence exists because u is upper- 
semicontinuous. We consider now the test function tp n — minl^,^™}, and notice 
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that (p n = ip in a neighborhood of (x, t), and hence u — tp n still has a maximum at 
(x,t). By the very definition of sub-solutions we have 

<p t (x,t) + F(x,D<p(x,t), D 2 cp{x, t),{y : <p n (y, t)>u(x,t) }) < 0. 

Consider K n = {(p n (-,t) > u(a;, t)} 2 -K" = > u(x, t)}. Since the sequence 

of the sets K n is nonincreasing, K n —¥ f\ .Kfc in the Kuratowski sense, and by 
construction K = f). Kk- It follows that 

liminf F (x, D<p{x, t),D 2 tp(x, t),{y : <p n (y, t) > u(x, t)}) 

n—>oo 

> F (x, Dip{x, t), D 2 <p(x, t),{y: u{y, t) > u(x, t)» 



and fl3.im follows. □ 

Remark 3.5. By the assumption iv) on F, a standard argument shows that if u 
is a subsolution (supersolution) and 9 : K — > R is increasing, then 8 o u is still a 
subsolution (supersolution) . 

3.3. A comparison result. Here we provide a comparison result, that is the main 
ingredient to get existence and uniqueness in the viscosity setting. Let us set 

Q T := R d x (0, T), d p Q T =R d x {0}, TZt — Qt U d p Q T - 

Moreover, we denote by USC(1Zt) and LSC(TZt) the space of upper and lower 
semicontinuous functions on IZt, respectively. The following comparison principle 
is an extension of [23j Theorem 1.7] for non local evolutions. 

Theorem 3.6. Let u € USC(1Zt) and v 6 LSC(1Zt) be a subsolution and a 
supersolution of (|3.8|), respectively. Assume that 



(3.12) limsup{u(z)-w(C) : {z, Q G d v Q T x K T ) U (K T x d p Q T ), \z- C| < r} < 0. 
Then u < v in TZt, and moreover, 

limsup{u(z) - v(C) : z,(e TZ T , \z - CI < r} < 0. 

rj.0 

Proof. The proof follows the line of the proof of [23, Theorem 1.7]. We do not 
provide a self-contained proof; we only indicate the changes needed to adapt that 
proof to the context of our non-local setting. For the reader's convenience, we will 
use the same notation as in [23], up to the fact that in our case, f2 = M. d (and the 
space dimension is denoted by d instead of N). 

As in [23] . by Remark 3.5 we may assume without loss of generality that u and 
v are bounded, and we extend their domain of definition on Q T by setting 

u(x, T) = lim sup{u(y, s) : (y, s) G TZ T , \y - x\ + \s -T\<r}, 

v(x,T) = lim ini{v(y, s) : (y, s) G TZ T , \y - x\ + \s - T\ < r} . 

The functions u and v are still upper and lower semicontinuous in Q Tl respectively. 

We first show that u is still a subsolution (and v a supersolution) in R d x (0, T] 
in the obvious sense. Assume indeed that u — ip has a strict maximum at z = (y, T), 
where tp is admissible and coercive. 
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Assume first that Dip(z) 7^ 0. For any n £N large enough, the function [x, t) M> 
u(x,t) — ip(x,t) — l/[n(T — t)] attains a maximum at a point z n = (y n ,t n ) £ Qt, 
where z n — ¥ z as n — > 00, moreover we have D(p(z n ) 7^ for n large. Hence, 

ft(z n ) + trr?_. s 2 + F ( X n, Dtp(z n ),D 2 (p(z n ), {(p{-,t n ) > <p(z n )}) < 0. 

Since any Kuratowski limit of {(f(-,t n ) > (p(z n )} is contained in {</?(•, t) > ip(z)}, 
using properties hi) and v.l) of F we deduce 

ipt(z) + F(x, Dtp(z), D 2 y{z), {<p(; t) > ip(z)}) < 0. 

If now D<p(z) = 0, we follow the lines of [23, Proposition 1.3]. Since (p is admissi- 
ble at z — (y, T), there are S > 0, / € T and u> € C°(R) with w(r) /r -> as r -> 
such that 

\<p(x, t) - <p(z) - <p t (z)(t - T)\ < f(\x - y\) + u(t - T) 

for all (x,t) £ B(z,S). Without loss of generality we assume that u £ C 1 (K) and 
w(0) = oj'(0) = and also that w(r) > for r / 0. Next choose a sequence 
oj n £ C 2 (M.) such that ui n (r) — > w(r) and oo' n (r) — > u>' n (r) locally uniformly in M. and 
set 

i>{x, t) = <p t (z)(t -T) + 2f(\x - y\) + 2u{t - T) , 

(3.13) $ n (x,t) = <p t (z)(t — T) + 2f(\x - y\) + 2uo n (t - T) - 1 . 

n(J — r) 

We have u — ip has a strict maximum at z. Hence for n large enough u ~ ip n has 
a strict maximum at z n — (y n ,t n ) £ Qt, with z n — > z, and tp n is admissible at z n . 
As u is a subsolution, we have, using also property iv) of F, 

(3.14) <p t (z)+2«/ n (t n -T)+ 1 

+ V \~f' F{y n ,y n - {M;t n ) > M^)}) < 
\Vn - y\ 

if y n ^ y, while tp t (z) + 2u' n (T - t n ) + l/[n(T - t n ) 2 ] < if y n = y. Letting n -> oo, 
we get tpt(z) < thanks to (|3.7[). Hence, as claimed, u is a subsolution in K d x (0, T}. 



Now, as in [23], we assume that 

(3.15) 6» := lim sup (u(z) - v(() : (z, C) G Q T , I* - CI < r\ > 

and try to get a contradiction. The proof then follows identically the proof in [23] 
from page 238 until the middle of page 241. In particular (using exactly the same 
notation), the case 9 — 6" is identical (since the non-locality does not play any 
role in that case), and we may jump to the case 6 < 6. As in [35], we then let 
(x, t, y, s) £ Q T x Q T be the maximum point of 

(3.16) u(x, t) - v{y, s) - af(\x - y\) - a(t - s) 2 -et-es- 5\x\ 2 - 5\y\ 2 , 

where / £ J, and e, a > are suitable positive constants, and 6 > is chosen 
in such a way that the value of this maximum point is strictly positive. We then 
can jump to the middle of page 241 (more precisely, up to "Now the definition of 
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viscosity solution yields"). Here, the situation is a bit changed. By Lemma 3.4 
get 

2a(i - s) + e + F fx, j|i + x + 2 ^I, {u{; i) > u(x, t)}j < 0, 

and 

2a(t -s)-e + F (y, af'Qp\)-^ - 26y, X - 251, {v(; s) > v(y, s)} \ > 0. 

Here X is a suitable symmetric matrix which also depends on S, and p :— x — y. 
Moreover, X, p, t and s are uniformly bounded, while \p\ is bounded from below. 
Therefore, we may assume that they converge, as 5 — > 0, to some limit denoted in 
|2"3"] by Y, p 0, i, s, respectively. Denote 

K & := {u(-,i) > u(x,i)} - x, L s := {v(-,s) > v(y,s)} - y. 

We may also assume that K$ — > K , L° s — » L c in the Kuratowski sense, for some 
K € C(R d ), L £ A(M. d ). We deduce (using the semicontinuity properties of F and 
the translational invariance) that 

2a(t-s) + e + F(o,af'(\p\)^- Y,k) < 0, 
(3.17) ) W { 

2a(i-s)-e + F{Q,af , (\p\)-^- Y,L) > 0. 



By (3.16) we have 

(3.18) u(x, t) - u(x, i) < v(y, s) - v(y, s)- 

(af(\x - y\) - af(\x - y\) + S\x\ 2 - S\x\ 2 + S\y\ 2 - %| 2 ) . 

Let R > and choose £ € K$ fl Br. Let also ij £ (0, 1/2) and q = 2r\p (recall 
p = x — y). Choose z with \z\ < i]\p\. Choosing x = x + £ and y = y + £, + q + z 
in ( 3.18| ), and observing that x — y = (1 — 2rj)p — z so that \x — y\ < (1 — we 
obtain since £ G ifa 

< w(i + £,*) — u(x, t) < v(y + £ + q + z,s) — s) — 
(a/(|p|) - a/((l - r/)|p|) - ^ ■ (24 + - J(£ + 9 + z) • (2y + £ + q + z)) . 

Since S(\x\ + \y\) ->■ (see [23]), |£| < i?, and 

af(\p\)-af((l-r,)\p\) >c> 

for some c independent of (5 (as \p\ is bounded away from zero), we have 



<*f(\p\) - - »?)|P|) - ^ • (24 + - + « + «) • (2$ + £ + q + z) > 

for 6 small. Thus, £ + q + z € La. As this is true for all |z| < 7j|p|, we find that for 
5 small enough, 

q+(K s nB R (0)) + B vm (0) C L A -. 

In other words, the sets g + (Lf«5 n L?r(0)) are at distance at least rj\p\ from L§. 
Taking the (Kuratowski) limits as S —> we deduce that dist(2ryp + K,L € ) > -q\p\, 
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and in particular that 2rip + K C L. Using property iii) of F and the translational 
invariance again, we deduce that 

F (-2 V p,af{\p\)^,Y,K) > F (o,af(\p\)^,Y,L 



\pv 'J- v ' J " v \p\" 

for any r\ £ (0, 1/2). Taking the limit r\ — > 0, and using the continuity of F with 
respect its first variable together with (3.17), we obtain that 2e < 0, a contradiction. 

□ 

3.4. Existence and uniqueness of viscosity solutions. To show the existence 
of viscosity solutions, we need the following stability result. 

Proposition 3.7. Let {u n ) n >i be a sequence of upper semicontinuous subsolutions 
of (3.8) and let, for any z — (x,t), 

u*(z) — limsup < u n (C) ■ \z — C| < T i n > - 
rio I r 



Then u* is also a subsolution of (3.8). 

Of course, a symmetric result holds for supersolutions. 

Proof. The proof of this result is a variant of the proof of [23l Prop. 1.3] (see also 
the proof of property (P2) in [H]), observing that if z n = (x n ,t n ) — > z — (x,t) 
and ip is a test function, then the sets K n := {ip(-,t n ) > (p(z n )} converge (up to a 
subsequence) in the Kuratowski sense to a set K C {f(-,t) > (fi(z)}. We conclude 
using the monotonicity and the semicontinuity properties of F. □ 

Given A C K", we denote by BUC(A) the space of bounded, uniformly contin- 
uous functions from A to ffi. We now can state a general existence and uniqueness 
result: 

Theorem 3.8. Let uq G BUC(WL d ). Then, there exists a unique viscosity solution 
u € BUC(M. d x [0,cx))) of (3.8) with initial condition uq. 

Proof. The proof of this result is very classical, see [TH1 [23] and based on Perron's 
method. We introduce 

u(x,t) = sup |u(x, t) : u subsolution of (3.8 ), minuo < u < maxtiQ , u(-, 0) < uo\, 

and its upper and lower semicontinuous envelopes. The fact that u* is a 

subsolution follows from Proposition 3.7 observing that at each point (x,t) we 
can find a suitable sequence of subsolutions (it„) n >i whose relaxed upper limit is 
u*(x,t). 

The fact that it* is also a supersolution is classical and obtained by contradiction, 
assuming that at some point z = (x,i) of (strict) contact with a test function tp < , 
ip does not satisfy (3.10). If Dip(z) ^ 0, one can use the test function <p to construct 
a new subsolution u > in a neighborhood of z, thus contradicting the maximality 
of u*. To treat the case Dtp(z) — one repeats the same construction, but (as in 
the proof of [23l Prop. 1.3]) with ip replaced by 

ip{x, t) = cp(z) + p t {z){t -t)- 2f(\x - x\) - 2uj(t - t). 



16 



A. CHAMBOLLE, M. MORINI, AND M. PONSIGLIONE 



The regularity properties of u and the fact that the initial condition is attained, 
can be shown as in the last part of the proof of [531 Theorem 1.8]. □ 

3.5. Application to our evolution problem. Here we show how to apply the 
viscosity approach developed above to our specific problem. First, we extend the 
Hamiltonian Ff in (3.2 1 to open sets, enforcing that the evolution of a closed set K 



agrees with the evolution of its complement, i.e., setting for every A 6 _4(R d ): 
(3.19) F f (x,p, X, A) := -F f (x, -p, -X, A c ) = F f (x,p, X, A c ), 



where the last identity follows by the very definition (3.2 1, (3.3) of Ff. It turns 



out, however, that the Hamiltonian Ff in (3.2 1 does not satisfy all the assumptions 



which are required in Theorem 3.8 In fact, it lacks assuptions v), v.l), v. 2). A 
basic counterexample is as follows: Let K be a, ball, and x n — > x 6 dK with x n £ K 
for all n. Then, F~ (x n ,p, X, K) = for all n and any p, X, while if p is the inner 
normal to dK, X is small enough and the radius of the ball K is large enough, 
then F~ (x,p,X, K) < 0. On the other hand, F+(x n ,p, X, K) will be constant 
(and positive). Hence Ff, in that case, will be l.s.c, but not continuous, and in 
particular v. 2) does not hold, neither the continuity with respect to x. 

In fact, we observe now that a continuous Hamiltonian extending the non-local 



curvature (2.181 does not exist. Indeed, let K = B := -Bi(O) C M 2 and x £ dB. Let 
moreover A n be open smooth subsets of B and x n € dA n be satisfying the following 
properties: 1) A n have vanishing diameter; 2) x n — > x; 3) the outer normal and 
the (euclidean) curvature of A n at x n agree with the outer normal and curvature 
of B at x, respectively. These conditions are clearly compatible. Set K n := B\A n . 
Then n^(K n ,x n ) — for n large enough (remember that /' = near 0). The idea 
now is that, if we could extend Kf into a semi-continuous Hamiltonian in the sense 
of v), it would follow that Kf(K,x) < 0, which is not true. More precisely, let 
u n = — djf n . By (3.4) and since Du n (x n ) = Du(x), D 2 u n (x n ) = D 2 u(x) we have 



(3.20) F f (x,Du(x),D 2 u(x),K) = Kt(K,x) > = liminf Kt(K n ,x n ) 

n— >oo 

= liminf Ff(x n , Du n (x n ), D 2 u n (x n ), K n ) = liminf F f (x n , Du(x), D 2 u(x), K n ). 

n— >oo n— >oo 

Since x n — ¥ x and K n — > K, we conclude that property v.l) does not hold. This 



means that Theorem 3.8 does not apply for our particular problem, without further 



smoothing (see Proposition 3.10 below). 

We can show a continuity slightly weaker than assumptions v), v.l), v. 2), which 
however will have some utility in the sequel. The following result shows that these 
properties are essentially true if {x,p, X, K) are of the form (x, D(p(x), D 2 ip(x), {ip > 
tp(x)}), when Dip{x) ^ 0. Finding a result similar to Theorem |3.8| but under this 



weaker assumption would be very interesting, and is a subject for future study. 

Lemma 3.9. Let (p n ,(p : M. d — > R be Cf oc functions, and assume that ip n — > tp in 
Cf oc as n — > oo. Let x £ M. d with D<p(x) ^ and consider a sequence (x n ) with 
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X n — > x. Then, 

(3.21) F f (x, D<p{x),D\{x), {if > ip(x)}) 



< liminf Ff(x n , Dip n (x n ), D 2 ip n (x n ), {(p„ > Pn(x n )}) 

and 

(3.22) F f (x, D<p{x),D\(x), {<p > <p(x)}) 

> limsup F f (x n , Dip n (x n ), D 2 ip n (x n ), {ip n > ip n {x n )}) . 

n—too 

Proof. We prove only the first inequality, the second one being a consequence of 
the first one and the identity 

F f (x, D<p(x),D 2 <p(x), > <p(x)}) = -F f (x, -D<p(x), -D 2 <p(x), {-<p > tp{x)}). 

First of all, replacing ip n with y H ► tp n (y — x + x n ) we may assume (by translation 
invariance of the Hamiltonian) that x n = x for all n > 1. Denote p = D(p(x), 
p n = Dcp n (x), p = Dip(x)/\Dip(x)\, and p n = Dtp n (x)/\Dip n (x)\, X = D 2 ip(x), 
X n = D 2 tp n (x), K — {ip> ip(x)}, K n = {ip n > ifn(x)}. One has that p„ -> p, etc, 
except for one detail: K n may not converge to K: more precisely any Kuratowski 
limit of a subsequence of K n is a set in between {if > p(x)} and K. 

Consider now s G [S, S] (recall that / is constant on [0, <$]). Since K is C 2 near 
x (by the implicit function theorem), there exists a positive s* £ {0,6], such that 
dist(a; — sp, K) = s if s < s* , and dist(a; — sp, K) < s if s G (s*, S] (possibly empty). 

We want to prove that for a.e. s (in fact, for all s ^ s*), 



(3.23) 

Since, clearly, 

\Pn\ 

2s 



Ks(K,x) < liminf kJ (K n , x) 



det 



\'P> 



2s 



det 



i- ¥] v P xv P 



we need to show (3.23) only when the right-hand side is zero, or more precisely, 
when dist(x — sp ni K n ) < s for infinitely many n > 1. In this case, let rik be a 
subsequence such that K„ k — > K (Kuratowski) and dist(a; — sp nk , K rik ) < s for 
all k. Let y nk € K nk such that \x — sp nk — y nk \ < s, and we can also assume that 
Vn k y S K C K. There are two situations: 

• either y ^ x, in which case, since \x — sp — y\ < s, we have s > s* . Since 
for s > s* k^(K,x) — 0, we conclude that (3.23) holds for all s ^ s*; 

• or y — x, in which case there exists € [x, y nk ] such that 

<Pn k ( x ) < <Pn k (yn k ) = fn k (x)+Dip nk (x)-(y nk ~x) + l-{D 2 Lp nk {z k ){y nk -x))-{y nk -x), 



hence 



(3.24) 



< Pn k ■ (Vn k -x) + -(X nk (z k )(y nk -x))- (y nk 



x). 



Now, 

,2 



> \X ~ Sp„ 



2.s 



I + s + T—\(Pn k ■ (Vn k ~X)), 
\Pn k I 
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(3.25) 



and hence, dividing by \y nk — x\ 2 we have 

2s (p nk ■ (y nk - x)) 



> 1 



\Pn k \ 



\Vn 



Set £ fc = (y nk - x)/\y nk - x\. Up to a subsequence £ fc -> £ _L p. By ( |3.24[ ) 
and (3.25) we conclude 



\P\ 



and in particular, I — s' /\p\VpXVp has a negative eigenvalue as soon as 
s' > s. It follows that s* < s and, again, we deduce (3.231 for all s*. 

Now, it remains to take the integral for s € an d it follows, using Fatou's 

lemma, that: 

F+(x,D<p(x),D 2 ip(x),{<p > ip(x)}) 

< limini Ft (x, Dip n (x), D 2 ip n (x), {(p n > <p n (x)}) . 

n— >oo 



In order to show (3.21 ), it remains to show a similar inequality for . 

This time, we need to show that for almost any s € [6, S] , if, for a subse- 
quence, Kj(K nk ,x) < 0, then Kj(K,x) must also take the value — l/(2s) det[J + 
(s/\p\)VpXVp] + . But it means precisely that dist(a: nfc + sp nk ,K^ k ) = s for all k, 
hence ip > ip(x nk ) on the ball of center x nk + sp nk and radius s. Passing to the 
limit, we deduce that tp > ^(a;) on the ball of center x + sp and radius s, so that 
dist(x + sp,K c ) > s. The thesis follows. □ 

Finally, we build an approximation of the Hamiltonian Ff which will fulfill the 
assumptions (i-v) of Section [3j To this purpose it is clearly enough to approximate 
F p for fixed p; a possibility is a follows, 



(3.26) F e (x,p,X,E) 



2p 



det 



\P\ 



VpXVp 



H £ (dist(x - pp,E) - p) 



\P[ 
2p 



det 



I 



\P\ 



VpXVp 



H £ {dist(x + pp, E c ) - p), 



where H e (t) is a continuous approximation of the Heavyside function, which is 1 
for t > 0, for t < — e, and nondecreasing. 

In that case, 
(3.27) 

d-l 



\F e (x,p, ±I,E)\ < 



\P\ 



IP 



P_ 

\P\ 



1 



\P\ 



d-lN 



< c(\ P \) ^ \ P \ 



2-d 



as \p\ — > 0, and in dimension d > 3 this Hamiltonian is indeed singular. The 
following result is straightforward: 

Proposition 3.10. The Hamiltonian F e satisfies all the properties required in 
Section [3l 



In particular, by Theorem |3.8| we deduce existence and uniqueness, in the vis- 
cosity setting, of the geometric flow corresponding to the regularized non local 
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curvature (3.26 1. In the next section we show an existence result (but with no proof 



of uniqueness) for the original non-local flow (3.51 



4. The geometric evolution associated to £f 
We will now follow a different approach, in order to construct a (level-set) flow 



of our curvature which is actually a viscosity solution of the equation (3.5), with Fj 



the Hamiltonian defined in (3.2 1. Here we assume that u Q is an initial datum with 
compact support, and bounded, uniformly continuous (uq € BUC c (R d )). 

The construction follows the approach first suggested by Luckhaus and Sturzen- 
hecker, and Almgren, Taylor and Wang [2"H H]. We follow here a simple strategy 
which has been elaborated in [TB] for the classical Mean Curvature Flow, and which 
we adapt to our setting. 

First, given a time-step h > and a compact set E, we define Tf~E (resp, T^E) 
as the minimal (resp., maximal) solution to 



(4.1) min I £ f (F) 



dist(x, dE) dx 



FAE 



£f(F) 



dE{x) dx 



dsix) dx, 



where c?_e(x) = dist(:r, E) — dist(x,]R ci \ E). The existence of a solution to (4.1) is 



not totally obvious, however, it can be established by considering the equivalent 
convex variational problem 

1 



mm < EHu, 

ueL 1 {M. d ;[0,X\) I h 



u{x)dE{x) dx 



with £f defined in (2.9), and observing that, given a solution of that problem, for 



a.e. s € (0, 1) the sets {u > s} and {u > s} are a solution to (4.1 ). The existence of 



a minimal (or maximal) solution follows from the fact that if E, E' are solutions, 



then also EnE' and E U E 1 are, thanks to (2.4). Moreover, it is not difficult to see 



that if F solves (4.1 ), then 



M f (F) =£ f (F), 



where A4*(F) is defined in (2.6). The following classical lemmas hold. 



Lemma 4.1. If E CC E' , then T+E C T^E' . Moreover, if E C E' , then T±E C 
T±E'. 

Proof. The proof is classical and we just sketch it. We first assume that E CC E' , 
so that ds > dE' a.e. We compare the energy (14. 1 h of F = E with the one of 



FnF', where F' = T~E', and the energy ([Ojl (with E replaced by E') of F' with 



the one of F U F' . We sum both inequalities and use (2.4) to deduce that F C F' 



Now, if ds > dE', we replace dE with + e and observe that the corresponding 
minimal solutions F E and F' satisfy F £ C F'. Let F be the Kuratowsky limit 
of F e (up to a subsequence). Then, it is easy to see that Fq is a solution, and 
T^E Cf Cf' = T^E'. The proof for T+ is almost identical. □ 

Lemma 4.2. Let E CC E' and let S = dist(dE,dE') > 0. Then T+E CC T~E' 
and, more precis ely, dist(dT+E,dT,;E') > 5. 
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Proof. Let z £ R d with \z\ < 5: z + E CC E' so that T+(z + E) C T^(E). By 
translation invariance of the scheme it follows that z + T^(E) C T7~{E), and we 
deduce the thesis. □ 

If E is a non-compact set with compact boundary, we can define T±E in a similar 
way (or simply let T±E = R d \ (T^(R d \ E))), and still the comparison holds. 



Thanks to the comparison lemma 4.1 starting from a function u £ BUC c (Mr) 
(with compact support, or constant outside of a compact set), for s > s' we have 
Tu {u > s} C T^~{w > s'}. It follows that we can define a function 

Thu(x) := sup{s : x £ 2^ {tt > s}} = sup{s : x £ T^{u > s}}. 

We easily see that for a.e. s, {T^u > s} = T^{u > s}. Using Lemma 4.2 we find 
that the distance between two such level sets of T^u is larger than the distance 
between the corresponding level sets of u: hence T^u £ BUC c (M. d ), with the same 
modulus of continuity. Finally, we can deduce (by approximation) that for any level 

s £ R, 

T^{u > s} = {T h u > s} , T+{u > s} = {T h u > s}. 
Now, starting from uq, we build a function Uh(x, t) : M. d x M. + —> R by letting 

u h (x,t) := (T h )^Uo 

where [•] is the integer part. By construction, has a uniform spatial modulus of 
continuity. The next lemma deals with the non-local evolution of balls. 

Lemma 4.3. Let x £ M. d , tq > and let Eg = B(x,ro). Then for every h, t > 
we have 

(T±)^(E ) = B(x,r±(t)), 

for some r^(t) > 0. Moreover, r^(t) r(t) uniformly in [0, T*(ro)] as h — > 0, 
where r is the solution to 

rS 

ds, 



(4.2) 



m= I f(s)[a + -) d - 1 -(a- -) + Y 

J 8 1 r r 

r(0) = r , 

andT*(ro) is extinction time of r(t) (i.e, such that r(T* '{tq)) = 0). Finally, there 
exists Co > such that for every ro < 1 we have 

(4.3) T*(r )>c r d . 

Proof. By translation invariance we may assume x — 0. Since the union of any 



family of minimizers of (4.1) is still a minimizer, we deduce that any rotation 
of (T^)^(Eo) is contained in (T£ )^ h \Eq) i.e., the maximal solution is radially 
symmetric. Analogously, by the stability of the minimality property with respect to 
intersection we deduce that (Tf^)^(Eo) is radially symmetric. By a rearrangement 
procedure it can be readily seen that the maximal and minimal solutions are in fact 
balls. Indeed, let r > be determined by \B(r)\ = (T+) lli] (E Q ). Then it is easy to 
see that 

£ f (B(r))<£ f ((T±)M(E )), 
/ dist(x, 8Eq) dx < / t dist(x,dEo) dx 
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with strict inequality whenever the radially symmetric set (Tj^)^ (Eq) is not a ball. 



For < r < R let e(r, R) be the total energy in (4.1 |for£ = B R and F = B{r) 
i.e., 



(4.4) e(r,R) = - j* f [s)w d [{r + s) d - ((r - s)+) d ] ds + ^ J\r - S ) S d " ] 



ds, 



where uJd denotes, as usual, the volume of the unit ball in M. d . A straightforward 
computation shows that J^e(r, R) — is equivalent to 

(4.5) Ur-R) = f S f'( S )[(l + -) d - 1 -((l--) + ) d - 1 ]ds. 

n Js r 1 

Now we construct the approximated evolution starting from B(r ). To this purpose, 
let us set rh,o = ?"o and define rh.% recursively, as the minimum point of e(r, rh,i-i) 
(and we stop if ru,% = 0). Denote by fh(t) the piecewise affine interpolation of r^h 
given by 

= r fc,[*] + (*~ ts]) ( r h,[i±i] - r M|]) ' 



Then, by ( [45) ), f h (t) solves 

d 



? M*) = *(Mm)); 

rh(0) = r , 



where 



5 



5 



(r):= / /'(«) (l + ^-i-^l- 



ds. 



r r 

Let [0,T*(r o )) be the maximal interval of definition for the solution to problem 



(4.2). Clearly, we have r(T*(ro)) = 0. Moreover, standards stability arguments 
in ODE yield that fh, and in turn r^, converge uniformly to r in [0, T] for every 
T < T*(tq). The uniform convergence in [0,T*(ro)] follows by monotonicity. 

Noticing that, for r < 1, |g(r-)| < cr 1_d for some c > 0, the final bound on T*(r ) 
follows by comparing with the solution to 

f r(f) = -cr 1 -^ 
\ r(0) = r . 

□ 

Lemma 4.4. There exists a time modulus of continuity w, such that for any 5 > 0, 
i/iere exists h(S) such that if x £ M. d , < h < /i(<5), and i, s > wit/i |i — s| < 6 
then 

\u h (x,t) - u h (x, s)| < 

Proof. Let wbea spatial modulus of continuity for uq, and therefore also for Uh(-,t) 
with i > 0. Fix r > 0. Then, Uh(y,t) < Uh(x,t) + w(r ) for all y £ B ra (x). 



Lemma 4.3 shows that if h is small enough, then Uh(x,t + s) < Uh(x,t) + U)(ro) 
for s < corg/2. Analogously, by Uh(y 7 t) > Uh(x,t) — o;(r ) for all y e B ro {x) we 
deduce Uh{x,t + s) > £) — w(ro) for s < co7"g/2. The thesis follows if we choose 

r = (25/co)V d ,uj(5) = w(r ). □ 
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Thanks to Lemma 4.4 we can extract a subsequence (hk)k>i such that Uh h 
converges locally uniformly in IR'' x M.+ to a function u(x, t) which is bounded and 
uniformly continuous in space and time. 

Remark 4.5. Let h n — > be such that Uh n admits a limit u. Then, as a straightfor- 
ward consequence of Lemma 4.3 we deduce that if for some level s£l, u(t, •) > s 
(resp., < s) on a ball of radius tq, then u(-, if) > s (resp., < s) on the concentric ball 
with radius r(t' - t), for t' > t, provided that r{t' - t) > (here r(-) solves (|4~2|). 

We can now show the main result of this section. 



Theorem 4.6. The limit u is a viscosity solution of (3.5), in the sense of Defini- 
tion 



Proof. First it is clear, by construction, that u(0, •) = uq. Hence we need to show 
that the equation holds for t > 0. We only prove that it is a subsolution, the 
proof that it is a supersolution being identical. Let ip g C°°(K <1 x R + ) and 
(x, i) € M. d x K + be a maximum point of u — <p. We may assume that this is a strict 
maximum point and that <p is coercive : if not, we should first replace (as usual) 
ip with ip(x, t) + r)(\x — x\ 2 + \t — 1| 2 ), derive an inequality for this modified function, 
and send r\ — > 0, which will give the desired inequality thanks to (3.21 1. 

By standard methods, we can then find (xk,tk) — > (x,t) such that tk > and 
u h k ~ has a maximum at 

Step 1. Let us first assume that Dip(x,t) ^ so that in particular, for k large 
enough, Dip(xk,tk) ^ 0. We have that for all (x,t), 

(4.6) u hk {x.t) < tp(x,t) +c k 

where c k := [u hk (x k , t k ) - ip(x k ,t k )], with equality if (x,t) = (x kl t k ). Let r] > 
and tp n hk : R d — > K given by 

= ¥>(tk,x) + c k + ^\x-x k \ 2 , 

then, for all x £ R d , 

u hk (t k ,x) < (fl k (x) 

with equality if and only if x — x k . Let e > and consider the open, nonempty 
set V e — {x : Uh k (t k ,x) > ip\ (x) — e}, which has positive measure, contains x k , 
and converges to {x k } in the Hausdorff sense as e — > 0. In particular, setting 
s s := Uf lk (x k ,t k ) — e/2, we have that for e > sufficiently small | W e | > 0, where 

VF e := u hk {t k ,x)>s e }\{x£R d : ^(x) > e + s e } CV e , 

Now, by minimality, we have 

£ f {{ u h k {-,t k ) > s e }) + t- I d {u h r.,t k -h k )>s e }(x)dx 

h k J{u hk (;t k )><>e} 

< £ { {{u hk (-, t k ) > s e } n {ifl k > e + s e }) 
1 



h k . 



d {u h A;t k -h k )>s c }{x)dx. 

{Uh k (-,t k )>Ss}n{!pl >S+S,} 
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Adding to both sides the term £f({v,h k (■, t k ) > s e } U > e + s 6 }) and using 

(|2~4l), we obtain 



Uh (■ ,t k -h k )>s s }( x ) dX < 0. 

'Ik JW F 



Observing that by ( [46] ), {u/ lfc (-, i fe - ft- fe ) > s £ } C {<£>(■> *k - fife) > s e - c fc }, we also 
have 
(4.7) 

d{tp{-,t k -h k )>s e -c k }[ x ) d>X < 0. 

Now notice that for z € W e we have 



(4.8) 

In particular, 
(4.9) 

Moreover, 



s e < (p(z,t k ) + c fe + -\z - x k \ 2 < E + s e 



W £ C B cv - e (x k ). 



(4.10) tp(z, t k - h k ) = <p(z, tk) - h k d t (fi{z, t k ) + hi / (1 - s)d 2 t tp(z, t k - sh k ) ds . 

Jo 

If y is the point closest to z with ip(y,t k — h k ) = s e ~ c k , so that \y — z\ = 
\d M ., tk -h k )>s e -c k }{z)\, then 

(4.11) tp(z, t k - hk) = <p(y, t k - h k ) + (z-y) ■ D<p(y, t k - h k ) 

+ f (l-s)(D 2 <p(y + s(z-y),t k -h k )(z-y))-(z-y)ds 
Jo 

= s e - cfe - d{ v (. ltk -h k )>8 S -c h }(z)\D<p(y,tk - hk)\ 
i 

(1 - s)(D 2 ip(y + s(z - y), t k - h k )(z - y)) ■ (z - y) ds . 

Combining pjjj ), ( |4.10[ ), and ( |4.11| ), we deduce 

} (z)\Dip(y,t k - h k )\ 

> -s + h k d t tp(z,t k ) - hi (1 - s)d 2 t <p(z,t k ~ sh k )ds 

Jo 

+ f {l-s)(D 2 ^(y + s(z-y),t k -h k )(z-y))-(z-y)ds. 



Note that, in view of (4.8), \ip(z,t k ) — ip(y,t k )\ <e + Ch k = 0(h k ), provided that 
e « h k are small enough. In turn, as \Dip{x k ,t k )\ ^ 0, we have \z — y\ — 0(h k ) 
and, using also (4.9 ), we deduce 

1 , , dt<p(z,t k )-% + Q(h k ) 

(4.12) Vk d w ,, tk _ hk)> _ s ^ Ck} {z) > ]D ^ tk ^ hk)l 

_ d t ip{x k , t k ) + O(Ve) + 0(h k ) 
~ \D<p(x k ,t k )\ + 0(h k ) ' 

We now focus on the term 

£ f {Wl k >e + s E }U W e ) - £f{Wl k > e + s E }) 
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of inequality (4.7 1. This is the sum of the two following expressions, which we will 
estimate separately: 

(4.13) / -f(s)(\(Wl > e + s E } UW E ) + B S \- |{^ > s + s E } + B s \) ds , 
Jo 



(4.14) / -na)(\{ V l k >e + 8 e }eB.\-\({ V l>e + a e }UW e )eB.\)da, 
Jo 

where A Q B denotes the set {x : x + B C A}. We recall that by assumption, 
f'(s) — for s < S, so that the integrals are in fact on [5, S). 

Let us first consider (4.13|. For any i in a neighborhood of x k , we have x G 
c?{t/?^ fc > <^ fc (x)} and we can define s*(x) G (0,5] such that for s G (0, s*(x)], 
dist(x + sv(x), {tp^ > <p1 k (x)}) = s and for s G (s*(x), S] (possibly empty), dist(a; + 
si/(x) 1 {<y5;[ fc > Lp^ k (x)}) < s. Here v(x) = —Dip 1 ^^ (x)/|Dc^ fc (x)\, and it is important 
to observe that thanks to the regularity of (p^ k , s* (x) is continuous near Xk ■ 

If e is small, for x 6 d{(p n hk > e + s £ }, there exists a minimal li (x) > 0, with 
h £ (x) < Cs/e, such that W E (1 {x + tu(x),t G [0,5]} C {x + tv(x),t G [0, 7i{x)]}, 
Clearly, for s > 6 and e small enough, 

(W v hk >e + s e }UW e + B s ) \ {{<pl u >s + s e } + B s ) 

3 {x + tv{x) : x G dWh k > £ + s E } ,s <t < min{s*(a;), s + h (x)}} 
The volume of this latter set is 

det(I + tVv(x))dtdH d ~ 1 (x) 



dWl k >e+s c }Jl(x) 

where I(x) is the interval (possibly empty) {t : s < t < min {s*(x),s + h (x)}. 
Fix <t > 0. A simple continuity argument yields that, for e sufficiently small, if 
5<s< s*(x k ) - ct, then s + \i (x) < s*(x). We deduce that 

(4.15) |(K fc > e + s e } U W e + B s ) \ ({^ >e + s £ } + B,)\ 

r- /•s + h'(x) 

> / det(J + tVv{x))dtdU d - 1 {x) 

J{p v h =e+s e , h s >0} Js 



S » det(I+(s + t)Vv(x)) J( _ M „, , 



det(I + sVis(x)) / ; . (T \ ' , , " dtdU d -\x) 

Wl k =e+s E ,h e >o} Jo det(/ + sVv(x)) 

> (det(J + aVv(x k )) + 0(Vi))(l + 0(>/i))|W e | , 

where we have used that 

\W e \ < / det(/ + tVi/(cc))dtd^ d - 1 (a;) . 

=e+s e , //>0} ./o 

Finally, integrating over s G [<5, s*(x fc ) — a] we deduce that (4.13) is larger than 

(4.16) (l + 0(Ve))\W e \J* iX } ff -(2 s /'( s ))( K +(^,a; fc ) + 0(^))d S 

= \W e \( K +(Klx k )+0(a)), 
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where the "curvature" (defined in (2.11) and (2.19)) is relative to the set KV = 
{Vfe ^ t fh k { x k)}- Here we have used the fact that y/s < a, so that 0(y/e) = 0(a) 



Now let us estimate the negative quantity (4.14). It is very similar, but not 
equivalent. 

We now need to introduce the function h E , defined for x £ ^{(p 7 ^ > e + s £ }, 
which is the largest real number (which can be nonzero only in a neighborhood of 
Xk) such that x + sv(x) £ W s for all s £ (0,h e (x)). For i in a neighborhood of 
Xk, we also introduce s*(x) £ (0, 5], such that dist(.T — sv(x), d{(p^ > (p 7 / lk (x)}) = s 
for s < s*(x) and dist(x — su(x),d{ip^ k > (p^ k (x)}) < s for s*(x) < s < 5. Just as 
a*, this quantity is continuous with respect to x, thanks to the regularity of <^ fc - If 
x £ d{ip^ > e + s £ } and s < s*(x), one has that 



u 



{<pl =e+se, h'>0} 



{x + tv(x) : -s< t < -s + h e (x)} 
2 (({< > e + s e } U W e ) e B.) \ (Wl h > e + s s } 9 B s ) 



at least when W e is small enough with respect to the scale s (which can be ensured, 
as we need to consider only s > 8). It follows that 

(4.17) | ((Wl > e + s e } U W e ) B s ) \ ({^ >e + s e } Q B s ) \ 

--s+h € (x) 

&et{I + tVv{x))dtd'H d - 1 (x) 



< 



{<p% =e+s s , h'>0} J- 



det(J - sW(x)) 

{ V l=e+s c ,h'>0} ./ii 



det(J- sVv(x)) 
< (det(J - sVv{x k )) + 0( v / i))(l + 0(y/e))\W t 



where we have used, this time, that 



\WA > 



/ det(J + tVv(x))dtdU d - 1 { 
Jo 



Fix it > as before. If s > s*(xk) + <J, then s > s*(x) + a/2 near Xk and we see 
that 

(4.18) | ((K fc > e + s s } U W £ ) G B.) \ ({^ > e + *J 9 S.) | = , 

provided that e is small enough. Integrating ( |4.17[ )-(4.18 ) over s £ {5,5} we deduce 



that (4.14) is larger than 



(4.19) (1 + 0(VE))\W S \(J -(2sf(s))( K j(Klx k ) + 0(V^))ds + 0(<7)) 

= \W e \( K J(Klx k ) + 0(a)). 

It therefore follows from ( |4.16[ ) and ( |4.19| ) that 

e f ({<pl>e + s e }UW e )-£t({<pl>e + 8 e }) > \W e \( Kf (Klx k ) + 0(a)). 
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Thanks to (4.12|, we deduce, after dividing (4.7) by \W S \ and sending e — > and 
cr — y 0, that 



(4.20) 



d t (p(x k ,tk) 



+ Kf(K%,x k ) + 0(h k ) < 



\Dip(xk,t k )\ 

where 0(hk) depends only on the regularity of ip. We may therefore send r\ to 
zero to deduce that (4.20) also holds also with KV replaced by K k = {ip(-,t k ) > 
<p{xk,tk)}i thanks to (3.21) and (3.4). Letting now k — > oo, using (3.21) and the 
monotonicity property iii), we deduce 

(4.21) d t <p(x,t) + F f (x,D<p(x, t),D 2 <p(x, *),{#,•) > <p(x,i)}) < 0, 
that is, u is a viscosity subsolution at (t, x). 

Step 2. Now we consider the case D<p(z) = and we show that <pt(z) < 0. Let -ipn 
be defined as in (3.13), with T replaced by t, and let z n — (x n ,t n ) be a sequence 
of maximizers of u — ijj n , such that x n — > x and t n — > ir . If x n ^ x for a (not 
relabeled) subsequence, then Dip n (x n , t n ) ^ and (4.21 ) holds for ip n at z n . Passing 
to the limit and using the properties of /, we deduce that <fit(z) < (see (3.14) for 
the details). 

We now assume that z n = (x, t n ) for all n sufficiently large. Set h n := t — t n and 



2h n 
Co 



where cq is the constant in (4.3). Note now that by ([3/7]) and (3.27), the function / 
appearing in the definition of ip n is of the form f(r) = g(r)r d , for a suitable function 
g such that g(r) — > + as r — > + . It easily follows that 



B(x,r n ) c {ip„(-,t„) < tp n (x,t n ) + f(r n )} 



2hr, 



Co 



C < u(-,t n ) < u(x, t n ) + g(r n ) 



2h n 
co 



Note that the last inclusion follows from the maximality of u — ip n at z n and the 
fact that u(z n ) = tjj n (z n ). By (4.3), the extinction time T*(r n ) of the ball B{x,r n ) 
under the non-local evolution satisfies T*{r n ) > c^r d = 2h n . Hence, by comparison 
(see Remark 4.5), we deduce that 



x e \ u(-,t) < u(x.t n ) + g(r n ) 



2K 

co 



Thus, using also the maximality of u — tp at z, 

tp(x, t n ) - ip(z) U(x,t n ) - u(x,i) 



h n h n 
Passing to the limit, we conclude that <pt(z) < 0. 



< g{r n )- 



□ 




Figure 2. A zebra and its smoothing (left: starting image, center: 
the non-local motion, right: the standard curvature flow), at a 
small time. 

5. Algorithm and numerical examples 

5.1. A numerical implementation of the time-discrete scheme. We show 
in this section an example of evolution with the motion studied in this paper, in 
dimension two. Actually, the implementation is not straightforward and only an 
approximate motion is computed, on a discrete rectangular grid. The approach we 
follow is described in |16j . It consists in minimizing, given a discretization of the 
signed distance function d E " 1 to the boundary of the set E 71 ^ 1 (negative inside, 
positive outside), the energy 

(5.1) minj(w) + ^rWu-d^W 2 

and define d E as the signed distance function to {u < 0}, computed as precisely 
as possible using a Fast-Marching algorithm [37J US]- Here, u, d E are defined on 
the discrete points ■ < i < N — 1 , < j < M — 1}, and the term 

h-d^f = <<':., r 

is the Euclidean norm. A spatial discretization term can be introduced in an obvious 



way. It turns out that if J is a correct approximation of the functional ( 2.2 1, then the 



algorithm is an approximation of the time-discrete scheme (4.1 ) studied in Section|4] 



In this case, the iterations should be an approximation of the motion driven by the 
energy. 

The discretization of the "total variation" J(u) is more complicated. Actually, 



the simplest here is to approximate (2.2) rather than (2.8). We fix p > 0, Let B be 
the discrete ball G Z 2 : i 2 + j 2 < p}, and let 



where the oscillation (here simply the max minus the min) is computed on the finite 
sets ((i,j) + B)n[0,N-l] x [0,M-1]. 

It turns out that for this particular energy, there is an approach, based on a graph 
representation and the maxflow/mincut duality, for minimizing binary problems 
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such as 



mm J{u) + y^fijUij 

i,3 



(given any real-valued matrix (fij)o<i<N, o<i<m)j and an algorithm for minimiz- 



ing (5.1) is easily derived. See [TB] for details and in particular [TBI Appendix B] 



for how this particular J can be implemented. 

5.2. Examples: two ways to shrink a Zebra. Figures [2j [3] and [2] show the 
motion applied to an initial set of curves with a lot of oscillations. As expected, the 
standard curvature motion shrinks the small scale objects much faster than the one 
based on the oscillation, in particular the stripes are preserved longer by the non- 
local flow. Notice that it is very difficult to estimate the exact corresponding times 
for the two flows, moreover, the numerical imprecision may provoke sometimes 
the "fusion" of the stripes in the classical curvature flow (wich is computed also 



using (5.1 1, but now J is a discretization of the standard total variation). 
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